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NARAYANA NUMBERS AND SCHUR-SZEGO COMPOSITION 



VLADIMIR P. ROSTOV AND BORIS Z. SHAPIRO 



Abstract. In the present paper we find a new interpretation of Narayana 

polynomials Nn{x) which are the generating polynomials for the Narayana 

- - - numbers N„ j. counting Dyck paths of length n and with exactly k peaks, 

OO ^ see e.g. 1161 . (These numbers appeared recently in a number of different 

^^ 1 combinatorial situations, [5l ll4|[T7] .) Strangely enough Narayana polynomials 

(_^ ' also occur as limits as n — ► oo of the sequences of eigenpolynomials of the 

Ch ' Schur-Szego composition map sending (n — l)-tuples of polynomials of the 

5_^ ^ form (x + l)"'~^(a: + a) to their Schur-Szego product, see below. As a corollary 

^ |. we obtain that every N„(x) has all roots real and non-positive. Additionally, 

.^r • we present an explicit formula for the density and the distribution function of 

' the asymptotic root-counting measure of the polynomial sequence {N„{x)}. 



1. Introduction 



(~| . 1.1. The Narayana numbers, triangle and polynomials. The Narayana num- 



bers N„^k, 1 < fc < n, apparently introduced by G. Kreweras in [llj are given by: 
where C* stands for the usual binomial coefficient, i.e. C* ''' 



?*■ I The latter formula immediately implies that for any fixed k the Narayana num- 

T~: ■ bers Nn^k are given by a polynomial in n of degree 2fc — 2 divisible by n. It is known 

^\ [ that Nn,k counts, in particular, the number of expressions containing n pairs of 

parentheses which are correctly matched and which contain exactly k distinct nest- 
ings and also the number of Dyck paths of length n with exactly k peaks. (Recall 
f^ ■ that a Dyck path is a staircase walk from (0, 0) to (n, n) that lies strictly above 

00 ! (but may touch) the diagonal y = x.) Some other combinatorial interpretations of 

^D ■ Nn,k can be found in [TBJ and references therein. 

The triangle 

rN ; 11 

c5' 13 1 (1) 

16 6 1 

1 10 20 10 1 

of Narayana numbers A^n.fc read by rows is called the Narayana triangle. (Later we 
will also interpret this triangle as an infinite lower-triangular matrix taking its left 
side of ones as the first column and its right side of ones as the main diagonal, see 

The generating functions of the rows of the above triangle are called the Narayana 
polynomials. More exactly, following the standard convention (see [2]) one defines 
the n-th Narayana polynomial by the formula 
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In what follows we will use the following notions. If P{x) is a univariate polyno- 
mial of degree n, then its reversion or the reverted polynomial P^{x) is defined as 
P (x) = x"P{l/x). A polynomial P{x) is called self-reciprocal if it coincides with 
its revertion up to a sign, i.e. P{x) — ±P^{x). Hence for any self-reciprocal P{x) 
if P{x) vanishes at xq then P{x) vanishes at l/a;o as well. A polynomial P{x) is 
called hyperbolic if all its roots are real. 

Remark 1. Each polynomial Nn{x) has a simple root at and each Nn{x)/x is 
self-reciprocal. 

The following simple 3-term recurrence relation satisfied by Narayana polyno- 
mials was found in [TBI P- 2]: 

(n + l)Nn{x) - (2n - 1)(1 + x)Nn-iix) - (n - 2){x - lfN„^2ix), (2) 
with the initial conditions Ni(x) = x, N2{x) = .x^ + x. 

1.2. Schur-Szego composition. The Schur-Szego composition (CSS) of two de- 
gree n polynomials P — J2^=oPj^'' ^^'-^ Q ~ S?=o IJ^'' ^^ defined by the formula: 

n 

When the same P and Q are considered as polynomials of degree n + k with van- 
ishing k leading coefficients then in accordance with the above formula one gets: 



Pn+kQ = Y.P^<l,^'l^n+f 



-1] 
3=0 



Extending these formulas one defines the composition of s polynomials by the for- 
mula: 

n 

Pi „;,. • • • „;, Ps - E^i.^- • • ■ps,j^'/ici+kr-' ■ 

3=0 

(For more details on CSS see [H [13].) 

Our main goal below will be a further study of a certain linear inhomogeneous 
map $„ initially considered in [SJ [Tj . Namely, in these papers the first author of 
the present paper has shown the possibility to present every monic polynomial of 
degree n with complex coefficients and vanishing at (—1) in the form: 

P = ^ai n ' ' ' n ^a„-i (3) 

where each composition factor Ka^ equals {x -|- l)"^^(a: + a^), a^ G C. (For the sake 
of convenience, we set K^o := {x + 1)"^^.) Now wc can introduce the map <&„. 

Notation 2. For any P{x) :— {x + l)(a;"^-'- +cix"^^ + • • • + Cn~2X + Cn-i) and for 
V = 1, ... ,n — 1 set Oy := X]i<ii<-<j„<n-i ^-h ' ' ''^j..' *-^- define a^, as the v-th 
elementary symmetric function of the roots of the composition factors presenting 
P{x). Finally, denote by <!>„ the mapping (ci, . . . , c„_i) i— > (ui, . . . , cr„_i). 

Obviously, <!>„ is linear inhomogeneous. The following theorem was proven in 

m- 

Theorem 3. (1) The mapping $„ has n—1 distinct real eigenvalues Ai^„ — 1, 

\ 71 \ n \ n^~ 
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(2) The corresponding eigenvectors are manic polynomials of degree n—\ van- 
ishing at (—1) and have the form: (x+l)"^^, x(x+l)"~^, a;(x+l)"'~^Qi^„(a;), 
..., x{x + l)Qn^3,n{x) whcrc degQj^nix) = j, j = 1,.. .,n-3, Qj>(-1) ^ 
0. The coefficients of each polynomial Qj,n{x) are rational numbers. 

(3) Each Qj.n{x) is self-reciprocal. More exactly, {Qj^n{x)) = {—^yQj.n{x). 

(4) The roots of each Qj,n{x), 1 < J < ft — 3, are positive and distinct. 

(5) For j odd (resp. for j even) one has Qj,n(l) = (resp. Qj,n(l) y^ Oj. 
Additionally, the middle coefficient in (a; + l)"~-'~^(3j^„(x) vanishes if n is 
even and j is odd. 

(6) For any j fixed and n —^ oo the sequence of polynomials Qj^n{x) converges 
coefficientwise to the monic polynomial Q*{x) of degree j which has ra- 
tional coefficients, all roots positive, and satisfies the equality {Q*{x))^^ = 
{—iyQ*{x) and the condition Q*(l) = forj odd. 

Remark 4. In Theorem\3[ we consider the action of $„ of the affine (n — 1)- 
dimensional space of all monic polynomials of degree n—1. If we extend this action 
to the ambient linear space of all polynomials of degree at most (n—1) then we 
acquire one more eigenvalue and eigenvector. Namely, the polynomial {x + 1)"^^ 
is the eigenvector of $„ ivith the eigenvalue 1 . 

1.3. Main results. Set Mj{x) — {—iy^^xQ*_i{—x). The most important result 
of the present paper (with a rather lengthy proof) is as follows, see details in 
Subsection [ 



Theorem 5. For any positive integer j the polynomial Mj{x) coincides with the 
Narayana polynomial Nj (x) . 

Part (6) of the above Theorem [3] then implies the following. 

Corollary 6. Narayana polynomials are hyperbolic for any n > 1. 

More information about the roots of Nn{x) is given below. For its proof consult 
Subsection Hm 

Theorem 7. (1) The number (—1) is a .simple root of Nn{x) for any positive 

even integer n. For n odd one has Nn{—1) ^ 0; 

(2) All roots of Nn{x) are distinct and nonpositive; 

(3) The roots of Nn-i{x)/x interlace with the ones of Nn{x)/x. Except for the 
origin the polynomials Nn-i{x) and Nn{x) have no root in common. 

Our final result is as follows, see details in Subsection 12.51 Given a polynomial 
P{x) of degree I define its root-counting measure jip ~ tX]j=i'^(-^ ~ ^0 where 
{xi, . . . ,xi} is the set of all roots of P{x) listed with possible repetitions (equal 
to the respective multiplicities) and 5{x — Xi) is the standard Dirac delta- function 
supported at Xi. Given a sequence {P„(a;)},degP„(a;) = n, n = 1, 2, . . . we call 
asymptotic root-counting measure of this sequence the weak limit /i = lim„^oo pLp^ 
(if it exists) understood in the sense of distribution theory. 

Theorem 8. The density p{x) and the distribution function k{x) of the asymptotic 

root-counting measure of the sequence {Nn{x)} of the Narayana polynomials are 

given by: 

11 2 

p{x) = — ; k(x) = 1 arctanV— a;, x <Q. (4) 

TT (1 — X)^/—X TT 

Remark 9. Notice that self-reciprocity of Nn{x) translates in the following (easily 
testable) property of p(x): 

.V«^.(i' 
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Figure 1. The theoretical distribution k(x) and the empirical dis- 
tribution of roots oi Niqq(x) on the interval [—1,0]. 
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2. Proofs 

2.1. Preliminaries. To prove Theorems [S] and [7] we will need a detailed study of 
the map $„ and, especially, of the equations defining its eigenvectors which in their 
turn give our polynomials Qj^n{x)- 

Notation 10. Set ek{j) = ^^(l, 2, . . . , j) to be the value of the k-th symmetric func- 
tion on the j -tuple of numbers (l,2,...,j), i.e. ek{j) = J2i<^i<...<^f,<j '^i ' ' ' ^k, 
k = I, . . ., j. Denote by cj)k{j) the sum 1*^ + 2*^ + • • • + j^ . 

Remark 11. The quantity efe(j) (resp. 4>k{j)) is a polynomial in j of degree 2k 
(resp. of degree k + 1) divisible by j{j -\- 1). 

Let Qj.n{x) :— x^ -f qix^^^ -!-•■■ + qj-ix + (—1)-' be the polynomial introduced 
in Theorem [3] and set Ij = (n — I) ■ ■ ■ {n — j) . Then by (1) of Theorem [3] one has 

Aj+2,n = n^ + V^j + 1- 

(the coefficients q^, depend also on j and n, but we prefer to avoid double indices.) 
By definition the polynomial Qj^n{x) satisfies the following relation: 

x{x + 1)"-^-2q^. „(a;) ^ 

= A,-,2,„x(a:+ir-i : {x+a,)[x+iT-' : ••• : {x+a,){x+iT~' I (x+ir-\ 

where {— ai, . . . , —Oj] is the set of all roots of Qj.n{x). After multiplication of both 
sides of the latter relation by Ij+i one gets that the coefficient i?^ of a;*^, fc > 1, in 
the right-hand side equals 

Rk := n^+'C^zl{C^,^,a, + C^zl) ■ ■ ■ {Ct^.a, + Ct-DC^-JiC^^y^' ■ 
The corresponding coefficient Lk in the left-hand side equals 



Lk := {n~l)---{n-j-l){{-iyct-U+Ct-]-2<lj-i+--- + C'^.-,-2<lj-k+i)- (5) 
Therefore one has q^ — Ov (see Notation [2]) and, finally. 
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(6) 
Thus the coefficients q^, i/— l,...,j— loi Qj^n{x) solve the system of equations 

(E) : {Lk = Rk, k = l,...,n-l}. (7) 

Lemma 12. The coefficients q^ can he expanded in convergent series: 

9. = ^) + ^ + r^ + --- (8) 

n — 1 (n— ly 

with respect to — ^j, where the numbers q^ € R are uniquely defined and indepen- 
dent of n . 

Proof: The coefficients q^ solve system (S). They are uniquely defined because the 
polynomials Qj^n{x) are uniquely defined by the eigenvectors of the mapping $„. 
They can be expanded in convergent series in -^—^ because the same property holds 

for the coefficients of system (S). As q^ are uniquely defined, thus q\f are also 
uniquely defined. □ 

Remark 13. We choose to expand q^, as a series in —^ (and not in —) because 
the eigenpolynomials of^n (see (2) of Theorem\^ are all of degree n — l. Besides, 
numerical computations show that it is the factor n—l and not n which appears 
most often in the denominators of the eigenvectors of the mapping $„ . 

Proposition 14. One has q[ = {—l)^qj2i = —j{j + l)/2. 

Proof: For k — 1 one has 

i-i 
Li = (-l)^(n-l)...(n-j-l), R^ = {n-l){l + Y.{n-lYq, + {-lY{n-iy). 

The equality Li ~ Ri can be written in the form 

(-!)■'■ (n - 1) • • • (n - J - 1) = i-iy{n - iy+^ + {n - ly qf2^ + o{{n - ly) . 

Hence qf\ + o(l) = {-lyUn - 1) ■ ■ ■ (n - j - I) - {n - iy+'^)/{n - ly . Observe 
that 

(n - 1) • • • (n - j - 1) = (n - iy+' - (1 + 2 + • • ■ + j)(n - 1)^' + o((n - 1)^) . 
The quantity qj2i depends on j, but not on n. Therefore {—iyq\_i = —(1 + 2 + 

■■■+j) = -j{3 + l)/2.a 

The next statement is central. 

Proposition 15. (1) For each (z^, ?) fixed the coefficient q{f is given by a real 

polynomial in j of degree 2{v + i). 
(2) For i = this polynomial is divisible by j(j + 1). 

2.2. Proof of Proposition [T5] . l". To prove part (1) of the proposition we use 
induction on i^ + i. Proposition 1141 constitutes the base of induction. The step of 
induction is explained in 2" - 3*^. 

Recall that the coefficients q^, give the unique solution to system (S). From 
now on we assume that system (S) is infinite, i.e. k = 1,2, . . .. Substituting the 
expansions ^ of the coefficients qi, in (E) we obtain a new system (denoted by 
(x)) with variables gi, , v = 1, . . . ,j — 1, i = 0,1,.... After this substitution 
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the equation Lk — Rk of system (S) transforms into an equation of the form 
Sl^; ^k,i/{n — 1)' = where the quantities A^^i are some Hnear inhomogeneous 

functions of the variables ql^ . (Notice that Ak.i depend on j but not on n.) The 
latter equation holds for all n e N if and only if all A^.i vanish. (The equation 
{Akj = 0} is denoted by {Ak,i).) 

2'^. The solution to system (x) is unique (which follows from the uniqueness of 
the polynomials Qj.ni^) for every fixed n, see Theorem |31). This solution depends 
only on j. The self-reciprocity of Qj^n{x) implies that q^ = (— l)^(7,y_^. 

In what follows we consider subsystems of system {x) of the form {{Ak.i),l = 
lo,. ■ ■ ,li}, i.e. systems defined in accordance with the filtration of the space of 
Laurent series in —^ by the degree of —^. We set I — s — j — k. 

n— 1 -J ^ n — 1 -J 

Notation 16. Denote by Ta,b the set of variables {qa , qa+ir ■ ■ v Qa+b}- 

To settle part (1) of Proposition[T5]we need Lemmas flTl and [TSl whose proofs are 
given after that of Proposition [121 

Lemma 17. The linear inhomogeneous form Ak,s-j-k depends only on the vari- 
ables in the set Jj-s,s-i := Ij-s,s-i ^Ij-s+i,s-2 U • • • UX,_i,o- 

3°. Suppose that the variables belonging to the set J'j-s+i,s-2 are already de- 
termined. (For s = 2 one has J'j-s+i^s-2 — ^j-s+i,s-2 — ^j-i,o ~ {QjJ-i}', see 
Proposition [HI) The system of s linear equations {A) := {{Ak,s-j-k),k — 1, . . . , s} 
is a system with s unknown variables, namely, those in the set 2j-s,s-i- This sys- 
tem has a unique solution (which follows from the existence and uniqueness of the 
polynomials Qj_n{x), see Theorem [2]). Hence the variables in the set lj-s,s-i are 
uniquely defined. 

Lemma 18. The solution to system (A) is an s-vector consisting of real polynomials 
in j of degree 2s. 

This concludes the proof of the step of induction in part (1) of Proposition [T51 
4^^. For V ^ \ part (2) of the proposition follows from Proposition [TH When 
solving the linear system [x) we express the variables in the set Tj-^.s-i as afiine 
functions of the ones in the set j7j_s+i.s-2- Suppose that all variables in that set 
are shown to be polynomials divisible by j{j + 1). Then the variables in the set 
2j_s,s-i will be divisible by j{j + 1) if and only if this is the case of the constant 
terms of system (A) (we call them CTs for short) . 

The CTs are the coefficients of [n — 1)"^^^'= of the expression (— l)^[/i — C/2 — 
{—\yU^ where 

t/i = (n - 1) . . . (n - J - l)Cl-_U ' f^2 = n^^'Cl-_\C':,_,{Clz\y , 

C/3 =n^"+iC^I^C„^i(C,ti)V(C,t)^+^ , 

see ([5]) and ([6]). In this difference the product U2 is irrelevant. Indeed, the highest 

power of {n — 1) multiplying any of the variables qi, in R^ is higher than the 

highest power of (n — 1) in C/2, see ©. 

Set (n-l)---(n-j-l) = {n-iy+'^ + V. Hence C/i = {{n-iy+'^ + V)Clz]^2- 
By Remark 111! the quantity y is a polynomial divisible by j{i + 1). Therefore for 
j = one has Ui = {n - l)C^Z\ = nC'^zlC^-i/Cl = C/3, and for j = -1 one has 
Ui — C^z\ = U3. Hence the CTs are divisible by j{j + 1). This completes the 
proof of Proposition [TS] Now we settle Lemmas [T7] and [T51 

Proof of Lemma [77[- 1°. Using the equalities C^zl/C^ = k/n and C^_i/C^ = 
{n — k)/n, one can present the equality Lk = Rk (see ([5]) and ([6])) in the form 

^n-l)...{n-j- m-iyctz]_2 + CtUl^-i + ■■■ + C°_,_2'7,-fc+i) = 
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= {n- fc)C^li(F + ^(n - kYk^-^q, + {-!)' {n - ky). (9) 

Replace in ^ the quantities Qj^ by their expansions ^. Consider the right-hand 
side Rk of ([9]) as a Laurent series in ^^ . Observe that if the integer k is bounded, 
then the foUowing relations hold: 

in-k)C^^Zl - ^"^,^^ +0{{n-l)'"^) , in-kyk^-'^ = fc-''-''((n-l)'^+0((n-l)'^-i). 

We use the last equality for v — j — s. The coefficient of (n — 1)^^''+'^ in Rf. is 
of the form: 

l(FgfJ, + k^-\fl,+, + ■■■ + kqf-,^'^ +H + r), 

where 7i is a linear form in the variables g/T with fi — m > j — s and r is a 
real number. Hence the form Ti contains only variables belonging to the union 
J7j_s+i_s_2 (because fJ. < j — I, m > 0) while the Hnear form {l/kl){k''qj2s + 

k'^~^qj-s+i'^' ' ' + f^1j-^i ) depends only on the variables in the set Xj_s,s-i- Hence 
Rk depends only on the variables in the set J^j-s.s-i- 

2°. Consider now the left-hand side Lk of ^. One can write 



B{n, j) := (n - 1) • • • (n - J - 1) = (n - iy+^ (l 



eiij) , 62 (j) 



see Notation [TUl For each i^ = j — k + I, . . . ,j the product B{n,j)C'^_^-'t2 ^ is a 



n — 1 (n — ly 



polynomial in the variable (ri — 1) of degree v + k. More precisely, 



(10) 

n-i-2 



-(g(''-^"+^)-q(''--'+^-i)ei(;.+fc-l)+...+(-l)''--'+^g(")e,_,+,(;.+fc-l)). 



Therefore, the coefficient of (n — \y ^^^ in the term i3(n, j)C^_-J^2 'J'l' of ^k is 
of the form 

1 

(v-3^k-\y: 

(11) 

The index v takes the values j — fc -|- 1, . . . , j — 1, see (0). Hence L^ is also a linear 
inhomogeneous form of the variables in the set Jj^s,s-\- ^ 

Proof of Lemma [TBI 1°. Consider equation (Ak^s-j-k)- Recall that its unknown 
variables are the ones in the set Xj-s,s-i- Present this equation in the form aiqj_g + 
■ ■ ■ + asqf_^ ^ (3 + G where the term /3 depends on j but not on the variables ql 

and t/ is a linear form in the variables gi- from the set J^j-s+i.s-2- 

2°. The quantity /3 is obtained by adding the terms (n — 1)^^ '*+'■" from the 
Laurent series of the expressions : A := {n — 1) ■ ■ ■ {n — j — l)(^l)"'C'^Ij-2' ^ -^ 



-(n - k)C'^zly and W := C,^Ii(-l)^(n - ky+'^ in equation ^. The coefficient 
of {n - 1).?-'*+'=^ in A equals (-l)^+^e,(j + k- l)/(fc - 1)! (see equation ^ with 
v = j) which is a degree 2s polynomial in j, see Lemma [TTl Its coefficient in W is 
a polynomial in j of degree s. Indeed, 

w = i-iyc^^zliin - 1) - (fc - i))^+i = {-lyc^zl £ ^Tf i(^ - i)"(fc - 1)'+'~^ 

7=0 
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IS 



where (— l)^C^_i is a polynomial in (n — 1) of degree k — 1 and Cj^'l'^ = C^^-^ 
a polynomial in j of degree s. As k < s < j, there is no term (n — lp^*+'^' in D. 
Hence /3 is a polynomial in j of degree 2s. 

3*^. The linear form Q is obtained from certain expressions in both sides of 
equality ([9]). First, one considers the terms -B(n, j)C^Zj_2~ Qi, in Lk and their 
coefficients of {n — 1)^^'+'^' given by formula (dJ). Recall that (see Lemma [T7|) 
the index i/ takes values > j ~ s. Set 6 := j ~ v. Hence 9 < s. In formula (fTTj) 
the term a := ql e^^j^^^s-s — ql e-s-e-s is the product of the degree 2(s — 9 — 5) 
polynomial eg-e-s in the variable j (see Remark [TT|) and of ±Qj_jy — ±<Ze which 
is a polynomial of degree 2{9 + (5) by inductive assumption. Thus a (and, hence, 
the whole contribution of L^ to the linear form Q) is a polynomial in j of degree 
2s. 

4°. Secondly, consider in R^ the product 



S:=C^Z\{n-kY+^k^ 



i^+ii.j-i'^ — 



9!^ 



GO 

Cl-_\k%{n - 1) + (1 - fc))''+i 5] g^") (n - 1)-" 

r;=0 



-,(')) 



Jri) 



The quantity qi, = ±qg is a polynomial in j of degree 2(0 + 77). 



Our goal now is to show that the coefficient of (n — 1) 



k+i^—r 



in the product 



C^_i{n — /c)""*"^ is a polynomial in j of degree r. (We prove this statement in 5*^ 
below.) This implies that the coefficient of {n — \y-~^+^ in 5 is a finite sum of 
polynomials in j of degrees r := 2{9 + rj) + r. To obtain a term [n — \y^^+^ we 
multiply the terms (n — 1)^+'^-'" and {n — 1)^''. In other words, one has {k + v — 
r) — rj — j- s + k, i.e. i' = j + r + ri — s and r = 2(j — i' + ri)+r = 2s — r< 2s. 
Thus the contribution of Rk to the linear form ^ is a polynomial in j of degree 
< 2s. The lemma is proved. 

5°. Proof of the latter statement. One has {n-kY+^ = ^^^^ C;+i(l - fc)''(n- 
^Y+i-r g^^^ Cl^i is a degree r polynomial in 1/, and therefore, also in the variable 
j—9 and thus in the variable j as well. The binomial coefficient C^Zi is a polynomial 

in (n-1) of degree (/c-1). Thus C^lJ(n-/fc)''+Ms of the form X;,^io '^r(n-l)'=+''-'' 
where dr is a polynomial in j of degree r. □ 

Now we finally return to our main results formulated in the introduction. 



2.3. Proof of Theorem [5j l". Consider the lower triangular matrix A4 whose j-th 
row contains the coefficients of the polynomial Mj{x) (starting with the coefficient 
of the linear term) followed by zeros. Let us turn the Narayana triangle ([T]) into an 
infinite lower triangular matrix of the form 



AA: 



/I 
110 
13 10 
16 6 10 
1 10 20 10 1 



(12) 



\: ■.:■.■. -.J 

Theorem El claims that the matrices A4 and Af coincide. Denote by Aii and A/j 
their l-th columns and by A4^ and A/"' their l-th diagonals (i.e. the sets of entries 
in positions {r,r + I — 1), r = 1,2,... in A^ and A/" respectively). 

The polynomials Mj{x) are monic and self- reciprocal with positive coefficients 
by definition. Therefore, one has Mi = Wi, M^ = N^ . 
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2". Suppose that the first m columns (and, hence, the first m diagonals as well) 
of A^ coincide with the first m columns (respectively, diagonals) of A/". The first m 
entries oi Mm+i and of Mm+i vanish. Their next m entries belong to the first m 
diagonals, hence, they coincide as well. 

3°. The entries of A^,„+i and the ones of A/'m+i (denoted respectively by Mm+i,j 
and Mm+i.j) are the values of polynomials R^^^ and i?™^^ in j of the same degree 
2m. For Aim+i.j this follows from Proposition [HI and for Mm+i.j it follows from 
the next formula for the Narayana numbers: 

^ Vm^m+i (j - l)(j - 2) ■ ■ ■ (j - m + 1) j{j - 1) ■ ■ ■ (j - m) 

] ■' m\ (m + 1)! 

For j = l,...,2m one has R^i^^U) = R-n'^^U), see 2°. For j = one has 



Nj,m+i = Mm+i,j = = Mm+i,j- The first two equalities follow from formula (flB]) . 
the last one can be deduced from part (2) of Proposition [T51 This proposition im- 
plies that Mm+i,j is divisible by j{j—l) (recall that Mj{x) — {—iy^^xQ*_-^^{—x)). 
Hence RZ'^\x) = R^^^x). D 

2.4. Proof of Theorem [7] . I*'. For n = 2 and 3 all statements of the theorem 
can be checked directly. Observe that Nn{x)/x has all coefficients positive. By 
Corollary [6] all roots of Nn{x)/x are real, hence negative, and is a simple root of 

Nnix). 

For any even n it follows from Nn{x)/x being self-reciprocal, of odd degree and 
with positive coefhcients that iV„(— 1) = 0. For n odd the polynomial Nn{x) does 
not vanish at (—1). This can be proved by induction on n. Namely, if Nn-2{x) 
does not vanish at (—1), then the same holds for Nn{x), see ([2]). 

2°. We prove the rest of the theorem using induction on n. Suppose that its 
statements hold for N2{x), . . ., Nn-i{x). Denote by ^i the i-th root of iV„_i(a;), 
£_i < Ci+i, ^n-i — 0. By part (3) of the theorem (proved for iV„_i(a;)) one has 
-^n-2(Ci)/6 7^ and the sign of Nn-2{^i)/£.i changes alternatively. By equal- 
ity ([2]) so does the sign of Nn{£,i)/£,i as well (equality ([2|) implies that the signs 
of Nn-2{S,i)/ii and iV„(Ci)/6 are opposite). As A^„-2(^«-2)/6i-2 > 0, one has 
Nn{^n-2) / £,n-2 < 0. The leading coefficient of the polynomial Nn{x)/x is positive, 
hence it has a root in (^„_2, 0) and by self-reciprocity a root in the interval (—00, ^1) 
as well. 

The polynomial Nn{x)/x is of degree {n — 1) and has at least one root in each of 
the intervals (—00,^1), (^1,^2), ■ • ■, (Cn-2,Cn-i)- Hence all these roots are simple 
(including (—1) for n even) and they interlace with the roots of Nn^i{x)/x. Thus 
the only root in common for Nn{x) and Nn-i{x) is which is simple for both of 
them, see part 1) of Remark [TJ D 

2.5. On root asymptotics of Narayana polynomials. In this subsection we 
prove TheoremEl Define *„(a;) = '^Ntix) ^^^ ^n{x) = „^"j(^j) , n = 1,2, . . .. Set 
^(a;) = liniji^oo "^nix) and Q{x) ~ lim„^oo 0n(a^) where these limits exist, ^(x) 
is called the asymptotic quotient and &{x) is called the asymptotic logarithmic 
derivative of the sequence {Nn{x)}. 

Lemma 19. The sequence {\f'„(a;)} of rational functions converges in C \ R<o to 
the function ^(x) = a; + 1 + 2y/x — (y/x + 1)^ where \fx is the usual branch of 
the square root which attains positive values for positive x. (Here R<o denotes the 
half-axis of non-positive numbers.) 

Proof of Lemma[TR- We need to invoke the classical result of H. Poincare, see [TJ p. 
287 - 298] and Theorem [^ and Remark ESl below. Indeed, dividing the recurrence 
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relation ([5]) by n + 1 we obtain the normalized reccurence 

2n — 1 n — 2 

^n(a^) -ri^ + ^)Nn-i{x) + -—-{x - l)^Nn^2{x) = 0. (14) 

n + 1 n + 1 

Taking limits of its coefficients when n — > oo we get from Poincare's theorem 
that the asymptotic quotient ^{x) for each x when ^(a:) exists should satisfy the 
following quadratic characteristic equation: 

*2(a;)-2(a; + l)*+(a;-l)2 = 0. (15) 

The exceptional set E d C (called the equimodular discriminant, see [3]) of those 
values of x for which the equation (|15p does not hold consists of all x for which two 
different solutions ^i(a;) and '^2{x) of (|15p have the same absolute value. 

Let us show now that in the considered case E — R<o- Indeed, two solutions of 
(fTS]) are given by ^^i (x) = a; + 1 + 2-^ and '^2{x) = a; + 1 — 2-^ for some choice of 
the branch of square root. One can easily check that if |^i(a;)| = |^2(a;)| for some 
value of X then a; + 1 is orthogonal to ^/x as two vectors in R^ ~ C. Denoting 
^/x = A + iB one gets x+l = A^ — B^ + l + i{2AB) and the latter orthogonality 
condition is given by the relation: 

A{A^ - B^ + 1) + B{2AB) = ^=^ A{A^ + B^ + 1) = ^=^ A = 0. 

But the real part A of ^/x vanishes if and only if a; is a negative real number 
implying E = R<o. 

Thus the asymptotic quotient '^{x) satisfies the equation (fT5| in C \ R<o. To 
show that \E'(a:) = x + 1 + 2^/x, i.e. it chooses the correct branch of solutions to 
([T^ we check that ^(1) — 1 + 1 + 2 = 4 and we prove that ^^(a;) is continuous in 
C \ R<o- 

Indeed, one knows that iV„(l) = X]fc=i ^n,k — Catn where Catn = ■:^nr[{^n) ^^ 
the n-th Catalan number, see 16J. Thus ^n(l) = n'^Ji) ~ ("+2)frt+i'i • There- 
fore, ^(1) — lim„^oo^n(l) = 4. In order to prove the required continuity (and 
analiticity) we use the well-known Montel's theorem claiming that a locally bounded 
family of analytic functions contains a subsequence converging to an analytic func- 
tion, see e.g [3j . For us it is technically easier to work with the sequence of inverses 
to ^n(ic), i.e. with the sequence {-jv-^^^t^}- We show that { at u\ \ is locally 
bounded in any compact domain separated from R<o. Thus since by Poincare's 
theorem { jv " "fi'i } ^^ pointwise converging in C \ R<o one gets by Montel's theo- 
rem that { jv {x\ } converges to an analytic function implying the same fact for the 
sequence {\I'„(a;)}. Indeed, by (2) and (3) of Theorem [7] all roots of each Nn{x)/x 
are simple, strictly negative and they interlace with that of A',i4-i(a;)/a;. Therefore, 
the partial fractional decomposition of ^ " %■) has the form 

n 



Nn{x) 

Nn+l{x) ^X-ai^n+1 



1=1 



where every ji is positive with X]i=i 7* — 1 ^^^d {ai_„+i, . . . , an,n+i} is the set of 
all non- vanishing (and therefore strictly negative) roots of Nn+i{x). If a; G C is 
any point denote by ^{x) its Eucledian distance to R<o- Let us check that for any 
a; £ C \ R<o and for any positive integer n one has 

Nnjx) ^ 1 

A^„+i(a;) ~ v{x)' 
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which immediately imphes the required local boundedness. Indeed, 



Nn{x) 



Nn+lix) 



v^ 7» 



•^ ^z,n+l 



< V ^' < Y- 7* 



■ ^ ^ ■T' — n.- .. , il f ^ 



i^l 



i,n-\-l 



I ^ iy{x) i/(x) ■ 



D 



Lemma 20. The sequence {9„(a;)} o/ rational functions converges in C\ R<o to 
the function Q{x) = —r-/=- 

Proof of Lemma \2(k Indeed, it is known, see |15j . that in case when both the 
asymptotic quotient ^{x) and the asymptotic ratio Q(x) exist and have continuous 
first derivatives in some open domain of C they satisfy there the relation 

ew - ^^ (16) 

A short sketch of its proof is as follows. Consider the difference 



^n{x) Pn+l{x) V "■/ Pn{x) 

Then one has. 



e(x) = lim -^ - e„(a;) = lim (n+1 



^(a;) n-»oo \^v]/^(x) / n^oo \(n + l)pn+i{x) npn{x) 

Applying p6p to the formula for '^{x) in Lemma [T^l we get 

D 

Notation 21. Given a finite measure fi supported on C define its Cauchy transform 
Cfj.{x) as 



C^(x) 



c a;-C 



The Cauchy transform of the measure is analytic outside its support, its -^- 
derivative coincides with the original measure and it has many more important 
properties, see e.g. [S]. Notice that for any polynomial P{x) of degree I the Cauchy 
transform C^p of its root-counting measure ^p is given by the formula C^p{x) = 
i.p(J\ ■ Therefore, given a polynomial sequence {P„(a;)}, degP„(x) = n,n = 1, 2, . . . 
we have that the Cauchy transform C^ of its asymptotic root-counting measure 
fi = Imin^oo fJ-Pn (if the latter measure exists) coincides with the limit C^(a;) = 
lim„^ ^^^^^ 



°° nP„{x)- 

The last result which we need to settle Theorem |5] and which is a particular case 
of Theorem 3.1.9 of [5] is as follows. 

Lemma 22. If pf^{x)^x € R<o is the density of a finite measure fi supported on 
R<o and Cf^{x),x G C \ R<o is its Cauchy transform then for any x G R<o one 
has 



Pt^{x) = ^ (Cm (a;) - Cp {xYj , 



where C^ (a;) = limt^a, C^ (x) with t belonging to the upper halfplane, C„ (x) = 
\mit-txCfi{x) with t belonging to the lower halfplane, and Z denotes the usual con- 
jugate of Z . 
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Finally, applying the latter formula to our case we get the required density 
formula in the statement of Theorem [HI Namely, 



P{x) = TT { — , ■ / -^= ) = Trr^ T ( (2; - iV^) -{x + iV^)) 

2tt \ X + i\/—x x — i\—xl l,ii(x^ — X) V / 



• —2i\/—x = 



27r(a;2 — x) " tt x^ — x 7r(l — x)y^—x 

Integrating the obtained formula for p{x) one gets the expression for the distribution 
function k{x) from the statement of Theorem [HI see Fig. 1. D 

3. Appendix. Poincare's theorem, 7, p. 287] 
1. Set-up. We consider a linear homogeneous difference equation of order k 

fit + k)+ aifit + fc - 1) + a2/(t + fc - 2) + • • • + akfit) = (17) 

with constant coefficients. Denote by 

Ai, . . . , A/c 

the roots of the characteristic equation 

A*" + fliA'^^i + a2\''^'^ -\ h flfe = 0, 

and assume that A^ have different absolute values denoted by ^i = |Ai|. We can 
then assume that 

6>6>--->6. (18) 

Since the roots Xi are all distinct then the general solution of (fT7|) is given by 

fit) = CiA* + C2A* + • • • + CkXl 

Let us choose one single solution of (fT7|) . i.e. assign some fixed values to the 
constants C^. Let Cp be the first nonvanishing among C^'s, i.e. 

Ci = C2 = • • • = Cp-i = and Cp ^ 0. 

Then one can show that 

t^oo fit) 

where fit) is the considered solution of p7p . Indeed, one has 

fit + 1) ^ CpA^+i + Cp+,\l\\ + ■■■ + CkXi+' 

fit) Cpyp + Cp+iXi^, + --- + CkXi 

For i — > cxD one has through (fT5|) that the limits of 

Ap+i \ / Ap+2 \ / Afc 



\ /'\ \ /''\X 
p / \ p / \ p 

equal 0, and we obtain 

hm l^tt^ = X 

t—>-oo 

that is the fohowing statement is vahd 



*— fit) 



Theorem 23. If fit) is an arbitrary (nontrivial) solution of the equation ( |J7| ) then 
the limit of the ratio f,^\ for t ^ 00 equals one the roots of the characteristic 
equation if only these roots have distinct absolute values. 
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We cannot say which root wiU be involved without the knowledge of the solution. 
One can only say that this root will have the number p of the first nonvanishing 
coefficient Cp in the considered solution. The following theorem of Poincare presents 
a generalization of this statement. 

2. Poincare's theorem. 

Theorem 24. // the coefficients Pi{t), i = 0, 1, . . . , fc — 1 of a linear homogeneous 
difference equation 

fit + k) + Pk-i{t)f{t + fc - 1) + Pk-2{t)f{t + k-2) + --- + Pa{t)f{t) = (19) 

have limits limj^oo Pi{t) — Oi, i = 0, 1, . . . , /c — 1 and if the roots of the character- 
istic equation 

A'^+flfc-iA'^-i + .-. + ao-O (20) 

have different absolute values then the limit of the ratio tif.) for t -^ oo of any 
solution f(t) of the equation il9\) equals one of the roots 

Ai, A2, . . . , Afe 

of the equation \20\) . i.e. 

Remark 25. In the present paper the role of the quantities f{t + i) is played by the 
values of the polynomials Nn for each x fixed. The parameter t is discrete - this is 
the index n. 
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